Experimental realization of directed percolation criticality in turbulent liquid crystals 
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This is a comprehensive report on the phase transition between two turbulent states of electro- 
convection in nematic liquid crystals, which was recently found by the authors to be in the directed 
percolation (DP) universality class [K. A. Takeuchi et al, Phys. Rev. Lett. 99, 234503 (2007)]. 
We further investigate both static and dynamic critical behavior of this phase transition, measuring 
a total of 12 critical exponents, 5 scaling functions, and 8 scaling relations, all in full agreement 
with those characterizing the DP class in 2+1 dimensions. Developing an experimental technique 
to create a seed of topological-defect turbulence by pulse laser, we confirm in particular the rapid- 
ity symmetry, which is a basic but nontrivial consequence of the field-theoretic approach to DP. 
This provides the first clear experimental realization of this outstanding, truly out-of-equilibrium 
universality class, dominating most phase transitions into an absorbing state. 



PACS numbers: 64.60.Ht, 64.70.mj, 47.27.Cn, 05.45.- 



I. INTRODUCTION 



Absorbing states, i.e., states which systems may fall 
into but never escape from, and phase transitions into 
them are expected to be ubiquitous in nature. For in- 
stance, spreading or contamination processes like fires or 
epidemics exhibit such transitions when the propagation 
rate changes: initially active (infected) regions eventu- 
ally disappear forever at low rates, i.e., the absorbing 
state is reached, while they can be sustained (pandemic 
regime) for fast-enough propagation. Examples abound 
far beyond: hundreds of numerical models, describing, 
e.g., catalytic reactions, granular flows, and calcium dy- 
namics in living cells, to name but a few, have been shown 
to exhibit such absorbing phase transitions p], Q . Such 
phase transitions also naturally arise from general prob- 
lems such as synchronization self-organized critical- 
ity H, spatiotemporal intermittency and depinning 
0, The vast majority of these transitions share the 
same critical behavior, that of the "directed percolation" 
(DP) class [3, [H, as long as they are continuous. 

Deep theoretical issues underpin this situation. 
Whereas universality is well understood for systems in 
thermodynamic equilibrium, this is still not the case for 
systems driven out of equilibrium, where even the rele- 
vant ingredients determining the class are not well un- 
derstood. In this context, absorbing phase transitions 
are central because of their genuine noncquilibrium char- 
acter, since absorbing states directly imply violation of 
the detailed balance. Janssen and Grassberger d E3] 
conjectured that the DP universality class contains all 
continuous transitions into a single effective absorbing 
state in the absence of any extra symmetry or conserva- 
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tion law. The DP class thus appears as the simplest, and 
most common case, as testified by overwhelming numer- 
ical evidence. 

However, the situation has been quite different in ex- 
periments. Over the last twenty years and more, a num- 
ber of experiments have been performed in situations 
where DP-class transitions would be theoretically ex- 
pected, but they have always yielded mixed a nd /or p ar- 
tial results with limited accuracy @, El El El El EE EE 
El El El M, M, HI (Table U). This lack of fully con- 
vincing experimental realizations, in contrast with the 
wealth of numerical results, has been found surprising 
and a matter of concern in the literature 0, [33] • 

Recently, though, studying a spatiotemporal intermit- 
tency regime occurring in turbulent liquid crystals, we 
found a transition whose complete set of static critical 
exponents match those of the DP class [3l[. The goal of 
the present paper is to provide a comprehensive report 
on this transition including not only a more complete de- 
scription of experiments designed to investigate the static 
critical behavior, but also new experiments giving access 
to dynamic critical behavior. 

The paper is organized as follows. We first illus- 
trate the coarse-grained dynamics of the turbulent regime 
of electroconvection studied and our basic experimen- 
tal setup (Sec. [TT|), with further details on image anal- 
ysis given in the Appendix. To characterize the critical 
behavior, we perform three series of experiments: (a) 
steady-state experiment under constant applied voltages 
(Sec. (b) critical-quench experiment starting from 

fully active initial conditions (Sec. IIV[) . and (c) critical- 
spreading experiment starting from a single active seed, 
prepared with a novel experimental technique using pulse 
laser, developed in this work (Sec.|V|. In Sec. IVI1 our re- 
sults are summarized and we discuss why clear DP-class 
critical behavior is observed rather easily in our system, 
contrary to many other experiments performed in the 
past in this context. 
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TABLE I: Summary of critical exponents measured in earlier experiments." 



(1+1)D System 


Size" 





v± 


t/M 




Mil 


Refs. 


annular Rayleigh-Benard 


22 




0.5 




1.9(1) 


1.9 


[11] 


annular Rayleigh-Benard 


35 




0.5 


0.5 


1.7(1) 


2.0(1) 


[12] 


linear Rayleigh-Benard 


25.7 


0.30(5) 


0.50(5) 


0.50(5) 


1.6(2) 


2.0(2) 


[12] 


interface roughening 






vx/v\\ = 


0.63(4) 






m 


viscous fingering 


60 


0.45(5) 






0.64(2) 


0.61(2) 


[13] 


vortices of fluid 


15 


0.5 








1.7 


[14] 


Taylor-Deen 


90 


1.30(26) 


0.64,0.53 0.73 


1.67(14) 


1.74(16) 


[15] 


Taylor-Couette 


70 


1 


0.4 




1.4-2.5 




[16] 


granular flow 






v \\ ~ Vl - 


= 1 






18, 23] 


torsional Couette 




u.ou^i ) 


0.53(5) 








flQl 
I19J 


rnvrnnninip cmVoc 
1CI I Oil muic opu^eb 


1 ns 




1.1(2) 


0.62(14) 1.70(5) 




1901 


lateral heat convection in annulus 


114 


0.27(3) 


0.30(4) 


0.75(3) 






[21] 


DP C 




0.276 


1.097 


1.734 


1.748" 


1.841" 


[24] 






v±/v\\ = 


0.633, 


V\\ — u± 


= 0.637 




(2+l)D System 


Size 





vx 


V\\ 




Mil 


Refs. 


liquid columns 


169 


0.56(5) 










[22] 


DSM1-DSM2 (present paper)" 


2.7 x 10 6 


0.59(4) 


0.75(6) 


1.29(11) 1.08(18) 


1.60(5) 




DP C 




0.583(3) 0.733(3) 


1.295(6) 1.204(2)" 


1.5495(10)" 


[25, 26} 



"Number in parentheses is the range of error given by the corresponding authors. 
"Number of effective degrees of freedom indicated by the corresponding authors. 

"Some exponents are obtained using scaling relations. Errors are then estimated by the law of propagation of error. 
"See also the remark [27ll . 

e Only some of the measured exponents are shown. See Table ITT1 for the complete list. 



II. EXPERIMENTAL SETUP 

We work on the electrohydrodynamic convection of ne- 
matic liquid crystals, which occurs when a thin layer of 
liquid crystal is subjected to an external voltage strong 
enough to trigger the Carr-Helfrich instability [H, l33j |. 
This is a priori a suitable system to study critical behav- 
ior thanks to its possibly large aspect ratio, fast response 
time, and easy controllability. We focus on the transition 
between two turbulent regimes, called dynamic scattering 
modes 1 and 2 (DSM1 and DSM2), observed successively 
upon increasing the root-mean-square amplitude of the 
voltage V at relatively low frequencies [H, [34[ . The dif- 
ference between DSM1 and DSM2 lies in their density 
of topological defects in the director field [Fig. [Ha)]. In 
the DSM2 state, a large quantity of these defects, called 
disclinations, are present [35[. They elongate and split 
constantly under the shear due to the fluctuating turbu- 
lent flow around. In DSM1, on the other hand, disclina- 
tions are present but kept smaller than the critical size 
and disappear immediately. Their density thus remains 
very low. The many disclinations in DSM2 lead to the 
loss of macroscopic nematic anisotropy and to a lower 
light transmittance than in DSM1. 

Our basic experimental setup is shown in Fig. [5] The 
sample cell is made of two parallel glass plates spaced 
by a polyester film of thickness d = 12 /im. Both in- 
ner surfaces are covered with transparent electrodes of 
size 14 mm x 14 mm, coated with polyvinyl alcohol and 
then rubbed in order that molecules are planarly aligned 
in the x direction, defined thereby. The cell is filled 



with A^-(4-methoxybenzylidene)-4-butylaniline (MBBA; 
purity > 99.5%, Tokyo Chemical Industry) doped with 
0.01 wt.% of tetra-n-butylammonium bromide. The tem- 
perature of the cell is kept constant carefully by a hand- 
made thermocontroller, composed of heating wires and 
Peltier elements controlled by a proportional-integral- 
derivative feedback loop with a lock-in amplifier [Fig. 
[Hb)] . Windows of the thermocontroller are made of sap- 
phire in order to improve the spatial homogeneity of the 
temperature. Throughout each series of experiments, the 
cell temperature is maintained at 26.0 °C with fluctua- 
tions typically of a few mK, unless otherwise stipulated, 
measured by three thermistors placed at different posi- 
tions close to the cell. 

We observe the electroconvection through the trans- 
mitted light from a handmade stabilized light source 
made of light-emitting diodes, recorded by a charge- 
coupled device camera. The observed region is a central 
rectangle of size 1217 /im x 911 /im (Fig. [21 inset). Since 
there is a minimum linear size of DSM2 domains, dj \[2 
[34j , we can roughly estimate the number of effective de- 
grees of freedom at 1650 x 1650 » 2.7 x 10 6 for the con- 
vection area and 143 x 107 ss 1.5 x 10 4 for the observation 
area. Note that the meaningful figure is that of the total 
system size, which is at least four orders of magnitude 
larger than in earlier experimental studies (Table [T| . In 
the following, we vary V and fix the frequency at 250 Hz, 
at roughly one third of the cutoff frequency 820 ± 70 Hz 
which separates the conducting and the dielectric regimes 
of electroconvection [H, [33| ■ 

Disclinations being topological defects, the sponta- 
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FIG. 1: (Color online) Spatiotemporal intermittency between DSM1 and DSM2. (a) Sketch of a DSM2 domain with many 
entangled disclinations, i.e., loops of singularities in orientations of liquid crystal. Blue dashed curves in the close-up indicate 
contour lines of equal alignment, (b) Snapshot taken at 35.153 V. Active (DSM2) patches appear darker than the absorbing 
DSM1 background. See also Movie SI in Ref. (c) Binarized image of (b). See also Movie S2 in Ref. (d) Sketch 

of the dynamics: DSM2 domains (gray) stochastically contaminate [c] neighboring DSM1 regions (white) and/or relax [r] 
into the DSM1 state, but do not nucleate spontaneously within DSM1 regions (DSM1 is absorbing), (e) Spatiotemporal 
binarized diagrams showing DSM2 regions for three voltages near the critical point. The diagrams are shown in the range of 
1206 /im x 899 )im (the whole observation area) in space and 6.6 s in time. 
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FIG. 2: (Color online) Schematic diagram of the experimental setup, in its entirety (a) and for the thermocontroller (b). LED: 
light-emitting diode, CCD: charge-coupled device camera, PC: computer, PID: proportional-integral-derivative. See text for 
details. 



neous nucleation of DSM2 in a DSM1 domain is in prin- 
ciple forbidden. It is indeed an essentially unobservable 
rare event, except along the edges of the electrode and for 
very high voltages, far from the range investigated in the 
present paper. Therefore, the fully DSM1 state serves as 
an absorbing state. On the other hand, DSM2 domains 
introduced externally, or present initially in the system 
can remain sustained in the bulk for large enough volt- 
ages, but eventually disappear for voltages lower than a 
certain threshold V c . Closely above V c , a regime of spa- 
tiotemporal intermittency (STI) is observed, with DSM2 
patches moving around on a DSM1 background [Fig. 



QJb), and Movie SI of Ref. 36]]. The basic dynamics of 
the observed STI is illustrated in Fig.QJd): active DSM2 
patches evolve in space-time essentially by contamina- 
tion of neighboring inactive (absorbing) DSM1 regions 
and by relaxation into the DSM1 state. This suggests 
an absorbing phase transition induced by change in rates 
of both elementary processes @ , which are functions of 
the applied voltage here. The order parameter p is then 
simply the ratio of the surface occupied by active DSM2 
regions to the whole area. 

Prior to any analysis, we must distinguish DSM2 do- 
mains from DSM1. This binary reduction can be easily 
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Voltage (V) 

FIG. 3: (Color online) The average DSM2 fraction p as a 
function of V in the steady state. Inset: same data in log- 
arithmic scales. Errorbars indicate the standard deviation 
of fluctuations in p(t) and V . Blue dashed lines are fitting 
curves. 



performed by our eyes, so we automated it, based on 
the facts that DSM2 domains have lower transmittance, 
longer time correlations, and have a minimum area of 
d 2 /2 [34[ (see Appendix for details). A typical result is 
shown in Fig. [TJc), and in Movie S2 of Ref. (36|. Figure 
[He) displays spatiotemporal diagrams obtained this way, 
showing how DSM2 patches evolve in the steady state. 
This supports the qualitative dynamics illustrated in Fig. 
[ljd) and indeed looks like the directed percolation of, say, 
water in a porous medium under gravitational field. 



III. STEADY-STATE EXPERIMENT 

We first observe STI in the steady state under constant 
voltage V, in the range of 34.858 V < V < 39.998 V. The 
voltage for the onset of steady roll convection (Williams 
domain) is V* — 8.95 V. Spatiotemporal distributions of 
DSM2 patches are recorded over the period 1000 s < T < 
8000 s, which is longer than 10 3 correlation times defined 
from the fluctuations of the order parameter pit). 

Figure [3] shows the time-averaged order parameter p, 
i.e., the average fraction of DSM2. It shows that the 
transition is continuous and that p scales algebraically 
near the critical point. Fitting these data with the usual 
scaling form [37J 



p~(v 2 -v c 2 ) 



2\0 



(1) 



the critical voltage V c and the critical exponent (3 are 
found to be 



V c = 34.856(4) V, = 0.59(4), 



(2) 



where the numbers in parentheses indicate the range of 
errors in the last digits [38(. Our estimate (3 = 0.59(4) is 
in good agreement with the value for (2 + l)-dimensional 
DP, /3 DP = 0.583(3) mill. 



We then measure the distributions of the sizes I and du- 
rations t of inactive (DSM1) regions, or intervals between 
two neighboring active (DSM2) patches. Histograms are 
made separately for each spatial direction, l x and l y , to 
take account of the anisotropy of DSM1 [3^, [4(|. Care 
is also taken to compensate missing intervals due to the 
finite observation window: since an interval of size I may 
not be captured within a frame of size L, i.e., either edge 
of the interval may not be in the frame, with probability 
l/L, the unbiased distributions Ni(l) are estimated from 
the observed ones by Ni(l) = N ohs (l)/(l - l/L) 0. 

The results are shown in Fig. [4] We find that the 
DSM1 distributions decay algebraically within the ob- 
served length/time scale for voltages very close to criti- 
cality, while for higher voltages they start to decay expo- 
nentially from certain characteristic length/time scales, 
as expected for continuous absorbing phase transitions 
0. The observed power-law decays are fitted as N\{1) ~ 
l-» x and Ni(t) ~ t~^\ with 

p x = 1.08(18), p y = 1.19(12), m = 1.60(5), (3) 

where p x and p y indicate the exponent p± measured in 
the x and y direction, respectively. These exponents are 
directly connected to the fractal dimensions df of the 
DSM2 clusters measured in the corresponding direction 
as p = d[ + 1. The critical exponents p± and p» for DP 
can be estimated using scaling relations with the order 
parameter exponent (3 and the correlation length/time 
exponents v±_ and v« , namely Q , 



fi± = 2 - @/v±, p\\=2- 0/v\\ 



(4) 



They give p^_ p = 1.204(2) and p^ p = 1.5495(10) for 
(2+1) dimensions [26|. Although the existence of slight 
discrepancies of order 10~ 2 in Eq. (H|) is suggested from 
numerical studies [H, [2i| , our estimates in Eq. ([3]) agree 
with values expected for the DP class at any rate. 

Moreover, the DSM1 distributions allow estimating 
correlation length and time scales, £j_ and respec- 
tively, from their exponential tails shown in Fig . JMd-f). 
Fitting distributions with an empirical form [111 [20j 



7Vi(Z),7Vi(t) ~ (Ar» + B)e 



-l/t 



(5) 



with powers p± and p» fixed at the estimates in Eq. ([3]), 
we obtain the results shown in Fig. [5[ Both £j_ and £y 
show algebraic divergence 

a~(V 2 -V c 2 )"^, £||~(V 2 -V C 2 )-^, (6) 

near criticality as expected, except for the first two points 
in Fig. [U^c), which deviate from the power law presum- 
ably due to finite length of movies (120 s) used to count 
the distributions. Fitting Eq. ([6]) to the data in scaling 
regions (the first five points for £ x , £ y and all the points 
except the first two for £||), we obtain 



v x - 0.75(6), u v = 0.78(9), 



1.29(11). (7) 
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FIG. 4: (Color online) Histograms of (inactive) DSM1 lengths l x ,l y and duration r in the steady state, in double-logarithmic 
scales (a-c) and in semi-logarithmic scales (d-f ) . The applied voltages are in ascending order from bottom right to top left (a-c) 
and from top right to bottom left (d-f), respectively. Dashed lines show the estimated algebraic decay at criticality. 
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FIG. 5: (Color online) Correlation length £ X) £ y and correlation time £y in the steady state, as functions of the deviation from 
criticality e = (V 2 — V 2 )/V 2 . Dashed lines are guides to eye. 



They are in good agreement with the DP values v± = 
0.733(3) and v^ p = 1.295(6) [Mill- In addition, the 
facts that no significant anisotropy is found between £ x 
and £y and that they can be even shorter than the cell 
depth d = 12 /zm [Fig. dja,b)] suggest that distributions 
of DSM2 patches are practically not influenced by the 
anisotropy of DSM1 [33, So| and by the existence of co- 



herent roll structure of width roughly d behind DSM1 

m- 

On the other hand, distributions of sizes of active 
(DSM2) patches do not show any long range correlations 
even close to criticality [Fig. EJa,b)]. It implies that the 
local recession of DSM2 into DSM1 sketched in Fig. QJd) 
is indeed always present. The effective relaxation rate can 
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/ x (^ m ) t (s) Voltage (V) 

FIG. 6: (Color online) Distributions of (active) DSM2 domain sizes in the steady state. (a,b) Histograms of DSM2 lengths l x 
in the x direction (a) and of DSM2 duration r (b). The applied voltages are in ascending order from bottom left to top right, 
(c) Characteristic length and time scales lxO,l y o,i~o, defined from the exponential decay in the histograms. The range of errors 
is smaller than the symbol size. Dashed line indicates the critical voltage V c . (Inset) Ratios of l x o and l y o to to (same symbols 
as in the main panel). 



be directly estimated, just as an inverse of the character- 
istic time To from the exponential tail of the temporal 
DSM2 distribution, shown in Fig. [6^c). The character- 
istic time To increases linearly with V except for high 
voltages, where it is not simply determined from the dy- 
namics of individual patches because of the saturation of 
the DSM2 fraction (Fig. [3]). In particular, it exhibits no 
sign of criticality. All these observations are consistent 
with, e.g., the dynamics of the so-called contact process 
[1^ (depicted in Fig. Did)), a prototypical model showing 
a DP-class transition [l|, 0] , indicating that such process 
indeed governs the coarse-grained dynamics of our STI 
regime in liquid crystal turbulence. 

Similarly, the spatial distributions of DSM2 sizes yield 
characteristic length scales I x q and l y o [Fig. [He)]. Al- 
though they appear to be also dominated by the re- 
laxation process, as suggested from an effectively con- 
stant ratio Io/to [inset of Fig. Etc)], they clearly show 
anisotropy [Fig. [5(c)], as opposed to the correlation 
lengths £ x and £ y estimated from DSMl-size distributions 
[Fig.[5j(a,b)]. This implies that the contamination process 
of DSM2 is indeed driven by the anisotropic, fluctuating 
shear flow of surrounding DSM1. The larger effective 
contamination rate in x is in line with the fact that the 
turbulent structure of DSM1 remains mainly in the x-z 
plane [39(, but at odds with the global elliptic shape of 
growing DSM2 nuclei, longer in y, observed for higher 
voltages [H, HH . Both would be explained if we assume 
that turbulent flow behind is faster and more correlated 
in the y direction, but further studies are necessary on 
this point. 



IV. CRITIC AL-QUENCH EXPERIMENT 

A typical experiment performed usually on numerical 
models showing absorbing phase transitions is the critical 
decay of active patches from fully active initial conditions 



[l|, |2j. In such critical-quench experiments, correlation 
length and time grow in time, and, as long as they remain 
much smaller than the system size, scaling estimates are 
free from finite-size effects. 

Experiments are performed as follows: we first ap- 
ply 60 V(^> V c ) to the cell and wait until it is en- 
tirely invaded by DSM2 domains. We then suddenly 
decrease the applied voltage to a value in the range of 
34.86 V < V < 35.16 V, i.e., near V c , and observe the 
time decay of activity for 900 s. We repeat this 10 times 
for each V and average the results over this ensemble. 

We first measure the decay of the order parameter p(t) 
after the quench [Fig. [7(a)]. As expected, p(t) decays 
exponentially with a certain characteristic time for lower 
voltages, converges to a finite value for higher voltages, 
and in-between at V = 35.04 V, it decays algebraically 
over the whole observation time. A simple scaling ansatz 
implies the following functional form for p(t) in this case: 

p{t)~t- a F p (et l l v V), a = P/v h (8) 

where e = (V 2 — V c 2 )/V c 2 is the deviation from criticality 
and F p (x) is a universal scaling function. From the slopes 
of the algebraic regimes for the three V values closest to 
the threshold, we estimate 

V c = 35.04(1) V, a = 0.48(5). (9) 

Note that V c measured here is slightly higher than in 
the steady-state experiments. In fact, the roll convec- 
tion onset V* — 8.96 V was also higher. We believe this 
is because of possible slight shift in the controlled tem- 
perature, and also of the ageing of our sample, a well 
known property of MBBA, during the days which sepa- 
rated the two sets of experiments. On the other hand, no 
measurable shift of V c was detected during a given set of 
experiments. We also confirmed that V c here is consistent 
with a threshold roughly estimated from steady state just 
before the critical-quench experiments. Our estimate of 
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FIG. 7: (Color online) Decay of the order parameter p(t) 
on critical quenching, (a) p(t) versus t, for V = 34.86 V, 
34.88 V, • • • , 35.16 V from bottom left to top right. The curve 
for V — 35.04 V (showing the longest power-law regime) is 
indicated by a thick line, (b) Same data with rescaled axes 
t[e|"H and p{t)t a , showing data collapsing. For V c , a, and v\\, 
we use values measured in the experiment [Eqs. Q and @] , 
but a collapse of similar quality is obtained also with DP- 
class exponent values. The dashed curve indicates the DP 
universal scaling function F® F (x) obtained numerically from 
the contact process. 



the critical exponent a is again in g ood agreement with 
the DP value cv DP = 0.4505(10) [if. 

Furthermore, the scaling form of Eq. © implies that 
the time series pit) for different voltages collapse on a 
single curve F p (x) when p(t)t a is plotted as a function of 
i|e|"ii. Our data do collapse reasonably well [Fig. 0(b)], 
where the upper and lower branches correspond to V > 
V c and V < V c , respectively. It is compared and found in 
good agreement with the universal scaling function F p (x) 
of the DP class [dashed curve in Fig. E|Jb)], calculated 
numerically from the (2 + l)-dimcnsional contact process. 
It shows that the decay of DSM2 patches is governed by 
this DP universal function, except for the very early stage 
where it is influenced by microscopic features of the liquid 
crystal, as expected. 




FIG. 8: (Color online) Autocorrelation function C(t,to) at 
V = V c , measured in the critical-quench experiments. Raw 
data in the inset are collapsed with rescaled axes t/to and 
C(t, £o)to' Dashed line shows the estimated asymptotic alge- 
braic decay. 

We also measure the autocorrelation function C(t,to), 
defined as 

C(t, t ) = (p(r, t)p(r, t )) - (p(r, t )) (p(r, t)>, (10) 

where (• • •) denotes average in space and over ensem- 
bles. During the critical decaying process, C(t, to) is not 
a function of t — to (inset of Fig. [8]) but scaled rather by 
t/to, a. feature sometimes referred to as "ageing" in the 
literature [44J. Our data can be collapsed (Fig. [5]) using 
the expected scaling form 

C(t,t )~to h F c (t/to), F c {x) ~ x~ Xc/z (z-kjo), 

(11) 

with 

6 = 0.9(1), X c /z = 2.5(3). (12) 

Both are in agreement with DP values 6 DP = 0.901(2) 
and A^ p /z DP = 2.583(14), estimated from the scaling 
relations (45| 

b = 2[3/v h \ c /z=l + {(3 + d s v ± )/v h (13) 

where d s = 2 is the spatial dimension. 

The dynamic aspect of critical behavior can also be 
characterized using first-passage quantities [4||. One of 
these is the local persistence probability P\(t), defined 
as the probability that the local state at a given point 
in space has not changed until time t, which typically 
shows a power-law decay P\(t) ~ t~ 01 . The non-trivial 
local persistence exponent 9\ is known to be in general 
independent of usual critical exponents such as (3, vx, uh, 
and far less is known about its universality, mainly due 
to the fact that persistence is a quantity involving an 
infinite-point correlation function. 

In the context of DP, local persistence is measured from 
the probability that initially inactive "sites" do not be- 
come active up to time t. (The persistence of activity is 
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mirror 



FIG. 9: (Color online) Local persistence probability P\(t) in 
the critical- quench experiments. The initial time is to — 2.14 s 
for the main panel and to = 0.81 s for the inset. The applied 
voltages are in ascending order from top to bottom. Dashed 
lines are guides to eye with the same slope. 



always dominated by the local relaxation into the inactive 
state and shows only exponential decay, similarly to Fig. 

Initial conditions are typ ically set to be random in 
numerical studies [47l IH, |49[ H(| , which is, however, im- 
possible in this experiment at present. We therefore con- 
sider, as "initial condition" for the local persistence, con- 
figurations chosen at some moment to during the critical 
quench and measure the persistence probability Pi(t— to). 

The scaling regime for the order parameter decay ob- 
served in Fig.[7]starting only after times of typically 1-2 s, 
we further limit ourselves to to values larger than this mi- 
croscopic, non-universal time. Typical results are shown 
for to — 2.14 s in Fig. [5J The local persistence proba- 
bility Pi(t) is found to converge to a finite constant for 
V < V c and to decay exponentially for V > V c , as ex- 
pected from numerical studies of absorbing phase tran- 
sitions 0, H El S At criticality, P x {t) decays 
algebraically (again beyond some time of the order of to) 
with 



01 = 1.55(7). 



(14) 



The value of 9\ for (2 + l)-dimensional DP is still a mat- 
ter of debate, with different estimates in past numerical 
studies: 6»f p = 1.50(1) ;48J, 6>p p > 1.6 and very re- 
cently 6f F = 1.611(7), found with an improved algorithm 
[Hlj . Our value, which is rather robust with respect to 
changing the initial conditions and the initial time (Fig. 

is in agreement with all numerical estimates to our 
accuracy. 



V. CRITICAL-SPREADING EXPERIMENT 

In order to complete the characterization of the dy- 
namic critical behavior of the DSM1-DSM2 transition, 
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FIG. 10: (Color online) Schematic diagram of the experimen- 
tal setup for the critical-spreading experiment. UV: ultravio- 
let, ND: neutral density. See text for details. 



we performed critical spreading experiments, which start 
from a single seed of active, DSM2 region. This allows 
to measure other critical exponents, such as the one gov- 
erning the scaling of the probability P s (oo) that a clus- 
ter starting from a single active seed survives forever: 
P s (oo) ~ e@ , known to serve as another order parameter 
characterizing absorbing phase transitions [l], Q . 

The two exponents [3 and /?' are known to be equal to 
each other for the DP class, thanks to an extra symmetry 
linked to time-reversal, the so-called rapidity symmetry 
[3, [|J . This symmetry also implies that a = 6 where 5 
is defined from the time decay of the survival probability 
at threshold, P s (t) ~ tr s . 

We stress, however, that the rapidity symmetry and 
the resulting scaling relations do not hold generically in 
absorbing phase transitions. Although a value of /3' con- 
sistent with DP was suggested from the critical scaling 
of hysteresis loops in the DSM1-DSM2 transition [H, 
it remains important to assess all spreading exponents 
as accurately and independently as possible, in order to 
complete our characterization of DP-class critical behav- 
ior and check directly the rapidity symmetry. 



A. Technique for nucleating DSM2 

In contrast to numerical simulations, it is not easy, in 
experiments, to prepare an initial, single, localized seed 
of DSM2 in an otherwise homogeneous DSM1 system. 
We developed an experimental technique for nucleating 
a DSM2 patch artificially, using a pulse laser. 

The experimental setup for the critical-spreading ex- 
periments is schematically shown in Fig. 1101 We emit 
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FIG. 11: (Color online) Probability of DSM2 nucleation in- 
duced by five successive pulses of laser. Laser is linearly 
polarized along or perpendicularly to the mean direction n 
of the molecules (squares and circles, respectively). Error 
bars denote the 95% confidence intervals assuming the bino- 
mial distribution. (Inset) Typical spatiotemporal diagram of 
nucleated DSM2 cluster for V = 36.65 V > V c . Scales are 
500 /im x 500 /im in space and 30 s in time. 

4-6 ns pulses of Nd:YAG laser (MiniLase II 20Hz, New 
Wave Research), focused by an objective lens (xlO, NA 
0.30, UPlanFLN, Olympus), into the cell. Using its 
third harmonic at 355 nm, around which MBBA has a 
broad absorption band [EH , l5~ij , and reducing its energy 
to roughly 0.3 nJ w 2 X 10 9 eV at the cell, we can indeed 
nucleate DSM2 from the absorbing DSM1 state (inset of 
Fig. ITTj) without any observable damage to the sample. 
Figure QT] shows the probability of the DSM2 nucleation 
induced by emitting five successive laser pulses at 20 Hz. 
The nucleation probability increases with voltage and is 
significantly higher for laser polarized along the mean 
director field. This confirms that nucleation is indeed 
brought about by the laser absorption of MBBA, since its 
absorbance is higher along the long axis of the molecule 
[53l [Hi]]. Moreover, the electronic structure of MBBA 
reveals that the ultraviolet absorption band stems from 
the C-N and C-C bonds between the aniline and benzyli- 
den rings and is strongly coupled with twist angles there 
[54l ]. We therefore speculate that absorption of ultra- 
violet laser pulses might lead to a sudden conformation 
change in the molecular structure, creating a topological 
defect and thus triggering a DSM2 nucleation. 



B. Results 

We perform the critical-spreading experiments with 
the technique above. For each measurement we emit ten 
successive pulses polarized along the molecules (in the x 
direction), with which DSM2 is always nucleated for volt- 



ages of interest. In order to make reliable statistics for 
the survival probability, it is necessary to repeat experi- 
ments at least hundreds of times for each voltage. This 
led us to improve further the temperature control of the 
cell, both in short and long time scales. This is achieved 
by placing the thermocontroller in a thermally-insulated 
chamber (Fig. H0|). made of plastic foam stage, wall, and 
ceiling, whose temperature inside is kept constant by cir- 
culating constant-temperature water. The cell tempera- 
ture is finally stabilized at 25.0 °C with fluctuations at 
most 2 mK over two weeks. This allows us to repeat the 
experiment 563-567 times for each voltage in the range 
of 36.25 V < V < 36.65 V. The roll convection onset 
and the critical voltage roughly measured in steady state 
were V* — 8.55 V and V c w 36.5 V, respectively. 

We measure not only the survival probability P s (t) but 
also the volume V(t) and the mean square radius R 2 (t) 
of DSM2 clusters, averaged over all the repetitions, even 
if the cluster dies before time t. The following relations 
are then expected from the scaling ansatz [l], Q : 

P s (t)^i- A F s (e^ii), 8 = p'/v h (15) 
V(t) ~t e F v (et 1/u »), = (d B v±-0-p')/v\\, (16) 
J2 2 (t)~t f F r (et 1 / , 'ii), C = 2/z = (17) 

where F s (x), F v (x), F T (x) are universal scaling functions. 

The experimental results are shown in Fig. [12J Ex- 
cept for R 2 (t) where all the data are almost overlapping, 
which is also typically the case in simulations of DP- 
class models, the data show opposite curvatures below 
and above a certain voltage. Seeking for the curve with 
the longest algebraic regime in V(t), which is statistically 
most reliable, we locate the critical voltage V c at 

V c = 36.45(2) V. (18) 

We then measure the three critical exponents <5, 9, £ from 
the algebraic regime for three voltages around V c , yield- 
ing 

5 = 0.46(5), 61 = 0.22(5), C = 1-15(9). (19) 

All of them are in good agreement with the DP expo- 
nents, S BP = 0.4505(10), 6» DP = 0.2295(10), and C DP = 
1.1325(10) [26]. This confirms the scaling relations ex- 
pressing the rapidity symmetry. It is further confirmed 
by plotting the ratio of the survival probability P s {t) 
to the DSM2 fraction p{t) in the critical-quench exper- 
iments [inset of Fig. [T27a)]. which shows that the two 
order parameters become asymptotically proportional to 
each other 

P s (t) ~ m 2 p(t) (20) 

with the coefficient m 2 w 3.2. 

We also tested data collapse [Fig. fl2T d-f)]. which pro- 
vides reasonable results given the limited statistical accu- 
racy of the data. The obtained scaling functions coincide 
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f|epl(s) t\e\ v U(s) t\E\ v W(s) 

FIG. 12: (Color online) Results of the critical-spreading experiments, (a-c) Survival probability P s (t), volume V(t), and mean 
square radius R 2 {t) of clusters started from a single DSM2 nucleus, for V = 36.25 V, 36.29 V, ■ ■ ■ ,36.65 V from bottom to top. 
The data for R 2 (t) are mostly overlapping. The curves for V = 36.45 V = V c are drawn with thick lines. [Inset of (a)] The ratio 
of the survival probability P a (t) from the critical-spreading experiments to the DSM2 fraction p(t) from the critical-quench 
experiments, both at criticality. The dashed line is a guide to eye, indicating roughly the asymptotic ratio of P s (t) / p(t) . (d-f) 
Same data with axes scaled after the expected scalings (|15[) -(|17 [ ). We use values of V c and critical exponents measured in the 
experiments. The dashed curves indicate the DP universal scaling functions F® p (x) , F® F (x) , F® p (x) obtained numerically 
from the contact process. 



satisfactorily with those of the (2 + l)-dimensional con- 
tact process [dashed curves in Fig. fT27d-f)]. As in the 
data collapse of p(t) [Fig. [T^b)], the collapsed data show 
typical time scales above which the corresponding quan- 
tities are governed by the DP scaling functions. 



VI. SUMMARY AND DISCUSSION 

We have performed three series of experiments, 
namely, steady-state experiments, critical-quench exper- 
iments, and critical-spreading experiments, to character- 
ize the critical behavior of the DSM1-DSM2 transition in 
liquid crystal turbulence. Table HT1 summarizes the main 
results. We have measured in total 12 critical exponents 
with reasonable accuracy, typically over a few orders of 
magnitude of power-law regimes. All of the measured 
exponent values agree within a few percent with those 
defining the DP universality class. Given that most of 
them are theoretically linked through scaling relations, 
we can equivalently say that we have experimentally con- 
firmed those scaling relations, 8 in total, that connect the 



measured exponents (Table Mil) . Among them, we have 
confirmed in particular the rapidity symmetry a = 8, 
providing also the asymptotic amplitude of the ratio be- 
tween the two order parameters, m 2 w 3.2. Moreover, we 
have also tested the expected scaling forms of Eqs. ©, 
(fTT]). (fl~5|) - (fr7|) through data collapse, and found them in 
good agreement with numerically-obtained DP universal 
scaling functions. Based on all these results, we conclude 
that the DSM1-DSM2 transition constitutes an unam- 
biguous experimental realization of an absorbing phase 
transition in the DP universality class. 

We now return to our initial remark concerning the 
surprising scarcity of experimental realizations of DP- 
class transitions (Table HJ. 

One central difficulty lies in the necessity to avoid, as 
much as possible, quenched disorder, which is known to 
be relevant. It is theoretically known that such disor- 
der does affect DP criticality and even destroys it for 
strong disorder [55]. Recent theoretical and numerical 
studies show that even weak disorder changes the asymp- 
totic critical behavior [5(| H3, HH, but the characteris- 
tic length/time scale separating DP and disordered, so- 
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TABLE II: Summary of the measured critical exponents. 



Exponent 




DSM1-DSM2" DPf25, 26] 


density order parameter 


P 


0.59(4) 


0.583(3) 


correlation length 




0.75(6) 0.78(9) 0.733(3) 


correlation time 


v \\ 


1.29(11) 


1.295(6) 


inactive interval in space 




1.08(18) 1.19(12) 1.204(2)' 


inactive interval in time 




1.60(5) 


1.5495(10)'' 


density decay 


a 


0.48(5) 


0.4505(10) 


local persistence 


Oi 


1.55(7) 


1.611(7) [51] 


ageing in autocorrelator 


b 


0.9(1) 


0.901(2) 




\c/z 


2.5(3) 


2.583(14) 


survival probability 


S 


0.46(5) 


0.4505(10) 


cluster volume 





0.22(5) 


0.2295(10) 


cluster mean sqr. radius 


C 


1.15(9) 


1.1325(10) 



"For ^tx an< l "J_ exponents measured in x and y direction are 
shown in this order. 
6 See also the remark 12711. 



TABLE III: Experimentally confirmed scaling relations. 



Scaling relations 



LHS RHS DPf25, 26] 



V±=2-l3/v ± a (in x) 1.08(18) 

(in y) 1.19(12) 

Ml, =2-/3/V 1.60(5) 

a = /3/i/|| 0.48(5) 

b = 2/3/^n 0.9(1) 

Xc/z = 1 + (J3 + dsv±)/v\\ 2.5(3) 

8 = /3/u }] 0.46(5) 

9 = (dsVx-20)/v\\ 0.22(5) 
C = 2/z = 2u±/vn 1.15(9) 



1.21(8) 

1.24(10) 

1.54(5) 

0.46(5) 

0.91(10) 

2.62(17)' 

0.46(5) 

0.25(11)' 

1.16(14) 6 



1.204(2) 

1.204(2) 

1.5495(10) 

0.4505(10) 

0.901(2) 

2.583(14) 

0.4505(10) 

0.2295(10) 

1.1325(10) 



"See also the remark [27l |. 
'The value of u x is used for /Ltj_ . 



called "activated" critical behavior grows fast with de- 
creasing strength of disorder [56|, H3| • It is therefore im- 
portant to work with systems made of macroscopic units, 
where quenched disorder is expected to be negligible. We 
consider that the quenched disorder in our system, which 
may take the form of inhomogcncitics in the electrodes 
or impurities in the sample, is also sufficiently weak. 

Still, our DSM1-DSM2 transition seems to be the only 
fully convincing one. Three factors explain, in our view, 
why our experiments provided such clear DP scaling laws. 

(a) Large system size and fast response. As already 
mentioned, one great advantage to work with electro- 
convection is that very large aspect ratios can easily be 
realized. The number of effective degrees of freedom of 
our cell is 2.7 x 10 6 , which is orders of magnitude larger 
than any earlier experiment (Table [J). This considerably 
suppresses finite-size effects and allows to observe scal- 
ing on several orders of magnitude. Similarly, the typical 
microscopic timescales of liquid crystals are very short 
(of the order of 10 ms), providing accurate statistics in 
reasonable laboratory time. 

(b) Almost perfectly absorbing state. The condition 
for being an absorbing state, i.e., that the system can 
never escape once it entered, appears to be somewhat 



too strict from the experimental point of view. Indeed, 
spontaneous nucleation of the active state seems to have 
be en p resent at least in some of the past experiments 
0, which blurs the critical behavior beyond some 
finite scales. 



On the other hand, our active state, DSM2, consists 
of topological defects, whose spontaneous formation is in 
principle forbidden. Of course this may occur in practice, 
as suggested from the observation of vanishing hysteresis 
at the DSM1-DSM2 transition but the nucleation 
rate remains so low that we cannot directly observe it, 
constituting an almost perfectly absorbing state. 



(c) Fluctuating absorbing state. In most earlier exper- 
iments and numerical studies, the absorbing state has 
been essentially fluctuation-free, or laminar. This is in- 
deed a natural choice suggested from the minimal the- 
ory of DP, and causes of course no problem in numerical 
studies. In experiments, however, such absorbing states 
may typically cause long-range effects through the rigid- 
ity of their laminar pattern and/or the propagation of 
soliton-like objects, which may further reduce the effec- 
tive system size and even break DP scaling 0, In 
contrast, our absorbing state, DSM1, is itself a fluctu- 
ating, turbulent state. Long-range interactions are then 
likely to be killed by the local turbulent fluctuations of 
DSM1, which may have led to the observation of clean 
DP critical behavior. It is noteworthy to remark that 
the Chate-Manneville coupled map lattice [5{|, a deter- 
ministic numerical model for absorbing phase transitions 
with a non-chaotic (laminar) absorbing state, exhibits 
non-DP critical behavior probably due to soliton-like ob- 
jects propagating through laminar regions. Interestingly, 
an elementary modification to render its absorbing state 
itself chaotic does lead to DP scaling [6(| ■ 



In conclusion, we have experimentally found that the 
DSM1-DSM2 transition in the electroconvection of ne- 
matic liquid crystals, showing an absorbing transition 
into spatiotemporal intermittency, clearly belongs to the 
DP universality class. Both static and dynamic criti- 
cal behavior has been investigated with the help of the 
newly developed experimental technique to create a seed 
of DSM2, confirming a total of 12 critical exponents, 5 
scaling functions, 8 scaling relations, and in particular 
the rapidity symmetry, in full agreement with those char- 
acterizing the DP class in 2+1 dimensions. We hope that 
this first clear and comprehensive experimental realiza- 
tion of a DP-class transition will trigger further studies of 
absorbing phase transitions and of related situations with 
genuinely non-equilibrium critical properties. In this re- 



spect the recent works of Corte et al. |61| and of Mangan 
et al. [62| are especially encouraging, concerning experi- 
ments and realistic situations, respectively, for absorbing 
transitions with a conserved field. 
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APPENDIX: BIN ARIZING IMAGES 

Every analysis presented in the paper is performed us- 
ing binarized images, where DSM2 domains are distin- 
guished from the absorbing DSM1 background. We give 
here a detailed description of the binarizing algorithm we 
used. 

The binarization is carried out in the following manner: 
(a) We prepare three successive images taken at 15 frames 
per second, and remove the inhomogeneity of the incident 
light intensity, (b) We then normalize the obtained inten- 
sity I of the three images with respect to mean (-Zdsmi) 
and standard deviation J/dsmi of the DSM1 intensity at 
a given voltage, namely I norm = (I - (Idsmi))/Wdsmi- 
Note that we can separately measure the intensity of 
the fully DSM1 state even above the threshold V c , since 
DSM1 always appears first when the voltage is applied, 
(c) Since DSM2 domains have lower transmittance than 
DSM1, we extract the regions where the normalized in- 
tensity is less than a certain threshold —Ith- Here we 
chose I t h — 1-5, determined so as to obtain a good 
agreement with direct visual observations, particularly 
in movies, (d) Taking into consideration that DSM2 do- 
mains move much slower than the local intensity fluctu- 
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